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The  variational  theory  of  three  dimensional  motion  of  curved  twisted  and  extensible  elastic 
rods  is  obtained  based  entirdy  on  the  kinematical  variables  of  position  and  rotations.  The 
constitutive  rdations  that  denne  the  resistive  couples  and  the  axial  force  as  gradients  of  the 
strain  energy  function  are  established.  A  candidate  for  the  strain  energy  function,  derived 
on  the  basis  of  dassical  assumptions,  is  presented. 


INTRODUCTION 


In  the  Historical  Introduction  of  the  "A  Treatise  on  the  Mathematical  Theory  of 
Elastidty,"  Love  (1892)  narrates  that  in  1742  Danid  Bernoulli  wrote  to  Euler  suuesting 
that  the  differential  equation  of  the  elastica  could  be  found  by  making  the  integral  of  the 
work  done  or  the  square  of  the  curvature  a  minimum.  Acting  on  this  suggestion  Euler  was 
able  to  obtain  the  differential  equation  of  the  dastica  and  the  various  forms  of  it.  Thus  the 
concept  of  the  strain  energy  was  bom  and  the  foundation  of  the  variational  theory  of  dastic 
rods  were  laid  out.  The  equilibrium  equations  that  were  much  later  devdoped  by  Love  are 
applicable  to  an  initially  bmt  and  twisted  rod. 

Our  aim  in  this  paper  is  to  establish  a  variational  formulation  for  the  title  problem  and  in 
the  process  infer  the  existence  of  the  strain  energy  function  and  determine  the  constitutive 
rdations  that  rdate  this  function  to  the  bending  and  twisting  couples  and  the  axial  force 
within  the  rod.  This  devdopment  together  with  equations  of  motion  and  the  geometry  of 
deformation  define  a  direct  appro^  and  an  exact  nonlinear  theory  for  the  three 
dimension^  motion  of  a  one  cumensional  dastic  medium  capable  of  resisting  bending 
twisting  and  extension.  Going  a  step  farther,  in  order  to  actually  constmct  an  explidt 
form  for  the  strain  energy  function,  we  enter  the  realm  of  hypothesis  and  use  Kirchhoffs 
description  of  deformations  in  a  tUn  rod.  This  view  enables  us  to  determine  a  strain 
energy  function  that  can  be  used  in  engineering  applications. 

The  recent  history  of  investigations  of  the  rod  theories  consists  of  devdopments  along  two 
separate  streams,  the  direct  approach  and  approximations  from  three-dimensional 
continuum.  In  the  direct  approach  a  one-dimensional  continuum  view  is  pursued  and  the 
medium  is  supposed  endowed,  in  addition  to  its  position,  with  vector  fidds,  the  directors, 
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that  are  to  be  interpreted  appropriately  to  define  bending,  twist  and  extension  properties  of 
a  rod.  This  approach  has  its  origin  in  the  work  of  E.  and  F.  Cusserat  (19091  and  numerous 
investigations  have  contributed  to  it  among  them  Naghdi  (1982),  Nagndi  and  Rubin 
(1984),  Whitman  and  DeSilva  (19701,  Green  and  Laws  (1966),  and  Eriksen  (1970). 
Extensive  investigation  into  the  qualitative  aspects  of  the  nonlinear  theory  such  as 
questions  of  existence  of  solutions  and  global  behavior  have  been  carried  out  by  Antman 
(1976).  His  basic  work  entitled  "The  Theory  of  Rods"  (1972)  describes  these  theories  both 
as  approximations  to  the  three-dimension^  continuum  theory  and  as  a  one-dimensional 
continuum  with  directors.  The  work  presented  here,  although  pertains  to  a 
one-dimensional  continuum  does  not  use  directors,  but  is  formulated  entirely  on  the  basis 
of  kinematical  quantities  consisting  of  the  position  vector  of  points  along  the  curve  of 
centroids  and  the  orientation  andes  of  the  cross  sections  of  the  rod  relative  to  a  fixed 
coordinate  system.  It  is  a  generalization  of  the  work  of  Tadjbakhsh  (1966)  in  which  the 
theory  of  planar  motion  of  the  extensible  elastica  was  described. 

The  history  of  construction  of  approximate  theories  in  the  context  of  three-dimensional 
nonlinear  continuum  theory  is  also  varied  and  to  it  many  investigators  including  some  of 
the  above  authors  have  contributed,  see  for  example,  Naghdi  and  Wenner  (1974). 


KINEMATICS 


An  elastica  is  a  nearly  uniform  slender  rod  of  finite  len^h.  In.  the  unstressed  state  the 
centroids  of  the  cross  section  form  a  space  curve  C  that  is  called  the  reference  curve  with 
an  arc  length  s.  The  orientation  of  the  principal  axes  of  the  aoss  section  vary  continuously 
along  the  rod.  This  means  that  in  the  unstrained  state  the  rod  has  arbitrary  twist  and 
curvatures.  With  respect  to  an  inertial  Cartesian  frame  x  the  position  of  a  point  s  in  the 
unstrained  state  is  denoted  by  X  =  Xi(s)ni,  i  =  1,2,3,  with  ni  being  the  dextral  unit 
vectors  of  the  frame  x. 


The  cross  sectional  area  can  be  slowly  varying  function  of  s  and  will  be  denoted  by  A(s). 
As  the  rod  deforms  the  curve  C  acquires  new  configuration  c  that  changes  with  time.  The 
arc  length  along  c  is  denoted  by  4  tbat  depends  upon  s  and  t,  i.e.  (  =  ^  (s,t).  The  position 
of  a  point  s  on  c  at  an  arbitrary  time  is  x  (8,t)  so  that 

x(8,to)  =  X(8),  ^(s.to)  *  8  (la,b) 

where  to  is  a  reference  time  at  which  the  rod  is  in  the  unstrained  state.  Also 


§  =  ('  =  (*'i  xl)* 


(2) 


where  prime  denotes  differentiation  with  respect  to  time  and  summation  over  a  repeated 
index  is  implied.  The  strain  e  is  defined  by 

e  =  - 1  (3) 

/ 

where  e  >  0  denotes  extension  and  e  <  0  contraction.  The  strict  positivity  of  ^  implies 
that  -1  <  e  <  a. 


3 


Attached  to  any  point  s  of  c  a  Cartesian  coordinate  frame  y  will  be  assumed  and 
will  be  refened  to  as  the  body  reference  frame.  The  coordinate  axes  of  the  body  frame  are 
Xii  Yv  Yi  Yi  axis  pointing  in  the  direction  of  increasing  s  and  y^  and  y,  being  the 

principal  axes  of  the  cross  section.  The  dextral  unit  vectors  of  the  y  frame  will  be  denoted 
by  e^,  i  =  1,2,3,  with  their  orientations  at  the  reference  time  t^  being  Ej. 

Denoting  by  lij(s,t)  and  Ljj(s)  the  elements  of  the  matrices  of  direction  cosines  of 
the  dextral  sets  e^  and  Ej  one  has 

n,  =  ly  ej  =  LyBj  (4) 

and 

®i  ~  I  ~  (S^jb) 

The  angles  represent  rotations  between  the  corresponding  pairs  of  E^  and  e^  when 

these  directions  are  assumed  to  issue  from  a  common  origin.  Fig.  1.  These  angles  are 
determined  through 

C08v>i  =  ®i  •  =*  Iji  Lji  « i=l|2,3  (sum  only  on  j)  (6) 

The  direction  cosines  1^  are  characterized  non— uniquely  by  three  orientations  angles 
and  ^3(s,t).  These  angles  can  be  selected  in  a  variety  of  ways  and 
represent  three  finite  rotations  about  unit  vectors  e^  or  n^.  If  these  rotations  are  properly 
selected  the  fixed  orientation  n^  may  be  brought  to  any  arbitrary  body  orientation  Oj. 
Kane  et  al.  (1983)  list  at  least  24  possibilities  for  the  order  of  rotations  of  the  angles  9^, 

9^  about  the  body  set  of  unit  vectors  or  the  fixed  set  of  unit  vectors  n^. 

Regardless  of  the  particular  choice  of  orientation  angles  the  angular  velocity 
u  —  0/^6^  of  a  cross  section^  element  Ads,  Fig.  1  of  the  rod  is  determined  uniquely  from 

=^igh  iig  ^ih  I  (*  =  If)  (7) 

where»7{j^  =  6^^  (c^^  +  l)/2,  (no  sum  on  i.j.k)  and  is  the  alternator  tensor  with  the 
non"*zero  components  —  ^331  ”  ^331  —  ^  1  and  ^^33  —  —  ^313  —  ~1* 

The  curvature  vector  K  =  Kiei  of  the  rod  can  be  defined  in  a  similar  way  with  K3 

r^resenting  twist  and  Kiand  Kirepresenting  bending  curvatures  about  the  principal 
directions  m  the  cross  section.  Using  the  dynamical  analog  of  E.I.  Routh,  Love  (1944)  has 
noted  that  if  the  firame  7  were  to  move  with  unit  speed  along  the  curve  c  such  that  at  any 
point  ^  of  c  it  has  the  orientation  of  the  7  frame  at  that  point  then  the  angular  velocities 
and  (J3  will  be  the  principal  curvatures  Ki  and  K3  of  the  rod.  Also  the  angular  velocity  U3 
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will  be  the  twist  curvature  K3  of  the  rod.  Thus  the  formulas  that  define  the  angular 
velocities  from  direction  cosines  can  be  used  to  determine  curvatures,  provided  time 
differentiation  is  replaced  by  differential  with  respect  to  Therefore  one  has 


where  differentiation  with  respect  to  ^  has  been  replaced  with  differentiation  with  respect 
to  s  and  curvature  parameters  ki  =  (l+e)Ki  is  also  introduced.  For  future  use  one  may 
note  the  formulas  for  derivatives  of  direction  cosines  lij  and  the  unit  vectors  ej. 

lij  =  (1+e)  fghj  iig  ~  (9^>b) 

/ 

ei  =  Ckji  kj  ek  ei  =  ckji  wj  ck  (10a,b) 

if  Ki  be  the  curvatures  of  the  rod  in  the  unstrained  state  (e  =  0)  then  firom  (8) 

/ 

Ki  =  T^igh  Lig  Lih  (11) 

Since  es  and  are  both  unit  tangents  to  the  central  line  one  has 

x'i  =  (l+e)li,  (12) 

We  assume  that  the  center  of  mass  of  the  aoss  sections  coincide  with  the  centroids. 
The  linear  and  the  centnd  angular  momentum  per  unit  length  are  then  given  by 

p=:pAxini  (13) 

and 

H  =  plu  =  p(Iii  oJiCi  +  Ijj  wj  Cj  +  Ijj  0/363)  (14) 

where  p  is  the  mass  density  per  unit  unstrained  length  and  I  is  the  diagonal  moment  of 
inertia  tensor  with  components 

I„  =  /y5dA,lM  =  /yidA,  =  +  (15) 


Referriim  to  the  body  set  of  axes  ei  one  can  define  the  vector  F  of  the  resultant 
shear  stresses  Fi  and  Fa  and  the  axial  stress  resultant  F3.  Similarly,  one  may  define  the 
couple  stress  vector  M  consisting  of  the  holding  moments  M|  and  Ma  and  the  torque  M3. 
Ex^dtly  we  have 


P  s  FiOi  and  M  s  Miei 


(16) 


The  wdl  known  dynamic  equilibrium  of  the  rod  can  be  expressed  by  the  equations 
of  the  balance  of  linear  momentum 
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F'+f=J  (17) 

and  of  the  balance  of  angular  momentum 

M'+x'»F  +  m  =  fi  (18) 

wherein  f  and  m  represent  distributed  force  and  moment  acting  on  the  rod.  The  scalar 
components  of  these  equations  can  be  referred  to  the  body  set  of  axes.  For  this  purpose  one 
needs  to  express  all  vector  quantities  in  terms  of  unit  vectors  e{  and  use  (9)-(10).  Then 
(17)  becomes 

f',  +  k^,-  k,F,+  =  pA  Iji  ('9a) 

f;  +  k,F,-  k,F,+  (5  =  M  ij  Ij .  (191>) 

Fi  +  k,F,-k/,+  t;  =  M*ilu  (19c) 

while  (18)  assumes  the  form 

M  j  +  kjM,  -  k,Mj  -  (l+e)Fj  +  m  ^  +  WjW,)  (20a) 

m',  +  k,Mi  -  kjMj  -  (l+e)Fi  +  mj  =  pl,(a;a  +  Wic;,)  (20b) 

M j  +  kjMj  -  kjMj  +  mj  »  pfJwj  +  (Ij  -  IJwjWj]  (20c) 

where  Ij  *  In  ,  Ij  *  I,,  and  I  =  Jjj  *  Ij  +  Ij.  The  superscript  y  on  the  components  of  f 
and  m  denote  the  components  of  these  vectors  in  the  body  reference  frame. 

To  express  the  equations  of  motion  in  the  inertial  frame  we  introduce  the 
components  of  the  stress  resultants  in  that  fnume.  Thus 

F*  s  ly  Fj  M]  *  ly  Mj  (21a, b) 

Then  (19)  becomes 

f]  +£^*pASi  (22) 

and  (20)  assumes  the  form 

M*  -(l+e)Fj  +  mj  *  />(I,j  ly  d;,  +  I,j  1^  w,w,)  (23a) 

Ma'-<l+e)F,  +  m ,  -  p(I,j  ly  w,  +  I,j  Ij,  w,a;,)  (23b) 
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M3  +  ni3  P(Irj  ^3j  *grj  ^sj  ^3f  ^s^r)  (23c) 

Either  of  the  set  of  equations  (19)Hf20)  or  (22)-(23)  can  be  considered  as  the  governing 
differential  equations  of  motion.  These  equations  will  have  to  be  supplemented  with 
constitutive  relations  that  define  resultant  axial  stress  Fj  and  the  resultant  bending  and 
twisting  couples  Mi,  M3,  Mj  in  terms  of  the  axial  strain  e  and  curvatures  k],  ks,  ks.  In  the 
next  se^ion  we  consider  the  derivation  of  these  constitutive  relations. 


I  —  I  2—2 

and  Pj ,  M{  and  P^ ,  are  the  applied  forces  and  moments  at  ends  s^  and  S3  respectively. 

-13  -12 

Also  V|’  and  are  the  initial  and  final  linear  and  angular  velocities.  The  strain  energy 


function  w  depends  upon  the  kinematical  variables  e  and  ki.  The  precise  nature  of  this 
dependence  is  the  constitutive  relations  that  we  seek  and  is  a  consequence  of  the 
stationarity  of  H.  The  functions  At  are  the  Lagrange  multipliers  that  allow  the  constraints 
(12)  to  be  incorporated  within  the  Hamiltonian.  As  a  result  xi  and  ly  can  be  regarded  as 
independent  variables.  Additionidly  the  constraint  (12)  implies  the  definition  j^-(3)  for 
the  strain  e  and  hence  in  (24)  e  can  also  be  viewed  as  an  independent  variable.  To  see  this 
we  need  to  note  that  if  eadi  side  cd  (12)  is  multiplied  by  itsdf  we  obtain  x'  x'^  = 

(l+e)^iilii  «  (l+e)>  which  is  restatement  of  (2)-(3).  The  terms  f^  Xj  and  m^  in  (25) 


represent  the  density  of  the  potential  of  the  applied  forces  and  moments  on  the  rod.  As 
stated  in  (5)  the  micpes  tp\  are  the  rotations  from  to  ei. 
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With  these  preUminaries  we  note  that  the  Euler  equation  corresponding  variations 
&ti  is  simply  =s  pA  which  when  compared  with  (22)  reveals  that  =  F^.  Next 

considering  the  variations  with  respect  to  e  we  obtain 

=  (2«) 

which  is  the  constitutive  relationship  determining  the  axial  force  Fj  as  the  derivative  of 
strain  energy  with  respect  to  axial  strain  e. 

We  now  turn  to  the  Euler  equation  corresponding  to  the  variation  5pi.  For  this 
purpose  we  note  that  lij,  ui  and  ki  depend  on  cpi.  For  orientation  angles  of  the  cross 
section  we  select  the  sequence  of  body  rotations  first  tfjea,  second  ^363  and  the  third  ^let 
with  0i  s  <pi.  In  this  sequence  the  last  rotation  is  through  with  respect  to  which 
variation  is  sought.  The  matrix  1  of  the  direction  cosines  is  given  by 

C3C3  — CjCaSs  +  SjSj  S1C3S3  +  CiSj 

1  =  B  ( W^3)A(tfi)  =  S3  CiCs  -S1C3  (27) 

SjC3  CjSaSs  +  SiCj  -S1S2S3+C1C2 

where 

Ci  =s  Cos  ,  Si  a  Sin  0{  (28) 

1  0  0  1  cos^  0  sin^ 

A(^  a  0  cos^  -sintf  B(tf)  a  0  1  0 

,0  sin^  cos^  -sintf  0  cos9 

'  ooiO  -9in9  O' 

C(^  a  sin$  009$  0 
0  0  1. 

Subsequently,  we  find  from  (7)  and  (8) 

S|  +  &!,  O^ja^]  C|C|  +  dj  Si,  ii/3=— SiC$  +  &sC| 
kfa^j  S$  +  ^1,  kaa^a  CtCj  +  tfi  Si,  kia— ^2  S1C3  +  ^3Ci 


(29) 

(30) 

(31) 
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From  (24)  we  have 


Noting  that  ffi  s  tpi  we  have  from  (30)  duifdipi  —  1,  duild<pi  =  uz ,  duzjdipx  = 

Also  from  (31)  dki(d(pi  =  0,  dkzfd<pi  =  kj,  dkz/dtfii  =  -kj  and  from  (27)  di\zld<px  =  -li2 
Using  these  results  and  the  inverse  of  (21a),  (32)  becomes 

”  ^*lk^a  -{l+e)Fa  +  m^  =  pli(^i  +  wjws)  (33) 

In  exactly  the  same  manner  one  can  proceed  o  determine  the  Euler  equation  for  a 
variation  Stp^.  Now  the  consecutive  sequence  of  body  rotations  ^jes,  9^\  and  ^262  is 
selected  with  9z  =  Without  going  into  details  one  obtains 

+  ni2  =  pl2(w2  +  U/iWs)  (34) 

For  variation  of  <pz  we  adopt  the  consecutive  sequence  of  body  rotations  ^lOi,  ^2^3 
with  9z  s  ^3.  The  matrix  of  direction  cosines  is 

C3CS  “C2SJ  S2 

l=sA(tfi)B(^3)C(tfs)=  S,S2Cj  +  CiSs  -SiS2S,+C,C8  -S1C2  (35) 

CjS2Cs  +  SjSs  C1S2SS+S1C3  C1C2. 
with  angular  velocities  of  the  cross  section  and  the  curvatures  given  by 

=s  ^iC2C3  +  ^283 ,  (1/2  *  —  ^^10283 ,  u/s—  ’9z+  &1S2  (36) 

ki  =  ^iC3C3  +  ^383 1  ka  *  9jCz  —  01C3S3 ,  k3  =  ^3  +  ^183  (37) 

For  this  case  lis  does  not  depend  on  ^3  and  hence  the  Euler  variational  equation  assumes 
the  form 


The  specified  boundary  conditions  at  s  s  S|,  83  must  be  consistent  with 


((Pi  -  f\  )Sxi  +  (liii  -  Mi)Ml!’  *  0  (39) 

■I 

for  arbitral^  and  independent  variations  ftci  and  6(pi.  Similar  restrictions  are  imposed  on 
initial  and  ual  data,  i.e. 
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A  STRAIN  ENERGY  FUNCTION 

In  oidet  to  gain  an  insight  into  the  nature  of  the  strain  energy  function  we  consider  the 
strain  of  the  lines  and  ai^es  in  the  cross  section  of  the  rod.  For  this  purpose  we  invoke  the 
Kirchhoff  hypothesis  which  assumes  that  the  plane  cross  sections  of  rod  that  are  normal  to 
the  axial  direction  in  the  unstrained  state  remain  normal  to  the  strained  axial  direction 
during  deformation.  Therdbre  the  position  vector  to  a  material  point  in  the  cross-section 
before  and  after  deformation  can  be  given  by 

R  *  X(s)  +  ytEt(s)  +  y,Ea(s)  (42) 

and 

r  =  x(s)  +  a(8)[y,ei(8)  +  yjCj(8)]  (43) 

respectively.  The  parameter  ofs)  is  to  be  fixed  by  enforcing  traction-free  boundary 
conditions  on  the  lateral  surface  ox  the  rod. 

Using  the  concept  of  extensional  strains  for  stretching  of  line  dements  and  distortion  of 
angles  between  perpendicular  lines  as  shear  strains  (Wempner,  1991),  we  define 
components  of  strain  by 

«ij  =  i(«i  •  81  “  Gi  •  Gj)  (44) 

where 

dt  *  * 

8a  *  ^  +  tt'yaea  +  oyiei  +  ogr jCa  +  (l+e)ea  (45) 

ffO  /  / 

G,  ■  “  yiBi  +  yaBa  + 


(46) 
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Using  (8)  we  can  establish 

Ci  *  Cj  =  Cijmka  (47) 

Bi  •  Ej  =  CijmKa  (48) 

where  Ki  is  the  curvatures  and  twist  in  the  unstrained  state.  Therefore  (9)  yields  as  the 
strain  components 

cii  =  f aa  =  ^(a*  —  1),  €ia  =  0 
ci3  =  i[a<a' yi  -  ayahs)  +  yjKs] 

C23  =  Wa'ya  +  oyjhs)  -yiKs] 

C33  =  e  +  ie*  -  yt[(l+e)aka  -  Ka]  +  ya((l+e)aki  -  KJ  -  yiy2(a2kik2  -  K1K2)  + 
iy?[a'»+a2(k*  +  k|)-K>-K>) 

+  \yl[a'  a  +  ai{k\  +  k|)  -K\-  K*)  (49) 

For  a  linear  isotropic  elastic  material  the  non>«ero  stress  components  per  unit  strained 
area  are 

<rit  *  (A  +  G)(a3  - 1)  +  Ac”  ,  cr”  =  0 

=  (A  +  G)(a3  - 1)  +  Ac** ,  a”  *  2G€as 

ffii  =  A  (a»  -  1)  +  (A  +  2G)«3« ,  a”  =  2Gc3i  (50) 

where  A  and  G  are  Lame’s  constant  and  the  shear  modulus,  respectively.  The  traction  per 
unit  und^ormed  are  is  then  given  by 

t3  =  aJigi  (51) 

One  can  define  the  adal  stress  resultant  Fs  by 

Fs  »  /*t*  *  63  dA  * 

•'1 

s  A  (1+e)  [(A+2G)(e+ie2)  +  A(o2- 1)] 

+  (A+2G)[l,I(l  +  e)ak,  -  Kjok,  +  laKl  +  e)Qka  -  Kajak, 

+  ^,(1  +  e)(a'  a  +  a^kj  -  KJ  +  aV  -  K|) 

+  ^3(1  +  e)(a'  a  +  aak*  -  KJ  +  oahj  -  K|)] 


(52) 
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Similarly  we  have 

Mj  =  J* oyjt*  •  ej  dA  =  a(A+2G)Ii^(l  +  e)aki  -  Kij  (1  +  e) 

+  IiQ!*ki^A(ar*  —  1)  +  (A  +  2G)(e  +  (j3) 

M,  =  —J  o^y^i^  •  Cj  dA  =  Of(A+2G)l2|^(l  +  e)ak2  -  K2j  (1  +  e) 

+  l2a»k2[A(a2  - 1)  +  (A  +  2G)(e  +  ie^)]  (54) 

Mj  =  -J* a(yjt®  •  Cj  -  yjt^  •  ei)dA  =  JGa*(a2kj  -  K3) 

+  Ja^ks  [a(  a2  -  1)  +  ( A  +  2G)(e  +  ie^)]  (55) 

One  may  note  that  the  integiability  conditions 


are  satisfied.  Hence  existence  of  a  strain  energy  function  is  assured  and  by  integration  we 
have 

W  =  A^^^^3(l+e+|^)  -  AA(l-aa)(e+ie2) 

+AIia2(o^l)^  +  (A+2G)liokip'^^|'^^Qkt— (l-i-e)Kij 

+Al20^fl*“l)^  +  (A+2G)l2oJr2^“^-^— ^~QJr7^1+e)K2j 
+^^l+e)a[l,(a'*  +  a^kJ-Kj)  +  Ha'i  +  aak!-Kl)] 

+AJa2(a2-l)^  +  ^^(e+ie2)(aak5-K!)  +  ^Ja(akr-K3)»  (57) 

For  an  initially  straight  rod  Ki  should  be  set  equal  to  zero.  The  above  form  of  W  reflects 
material  isotropy,  i.e.  W(e,ki,  k3,  ks)  =  W(e,  k2,  ki,  ks),  provided  that  Ii  =  I2. 
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We  note  that  positive  curvatures  imply  positive  bending  moments  and  conversely  negative 

curvatures  imply  negative  bending  moments  provided  that  e>(-l+l/fJ)  for  a  =  1.  This 
shows  that  equations  (53)  have  a  limited  range  of  validity  if  the  sense  correspondence 
between  moments  and  curvatures  is  to  be  retained.  We  also  note  the  second  order  coupling 
between  the  squares  of  the  curvatures  and  the  axial  strain  e  in  (52),  which  implies  that 
axial  force  can  be  generated  by  bending  or  twist  only. 

The  parameter  a(s)  depends  on  the  boundary  conditions  applied  at  the  lateral  surface  of 
the  rod.  If  the  lateral  surface  is  fixed,  then  afs)  =  1.  For  zero  tractions  on  the  lateral 
surface,  a  condition  appropriate  for  thin  flexible  rods  is  adopted  according  to  which  the 
average  of  cru  and  <7^2  over  the  cross-section  should  vanish,  i.e. 


J^aUdA  =  j[a"dA  =  0 

The  above  condition  reduces  to 

H(a,ki,e)  =  a'J  +  a2(Iiki  +  hk]  +  3kl  +  ^) 


(58) 


— (IiKt  +  IjKj  +  JK3)  +  2Ae(l  +  je)  =  0  (59) 

This  equation  should  be  interpreted  as  a  differential  equation  for  a(s)  when  ki  and  e  are 
known.  To  achieve  this  (52)  is  solved  in  an  iterative  procedure  in  which  at  every  step  ki 

and  e  are  known.  We  begin  by  writing  (52)  as  Um  H(aQ*j,  k^,  ej  =  0  with  n  =  0,  1,  2,  3, 

...  For  the  first  iteration  (nsO),  eo  =  0,  kf  s;  K^,  s  1,  and  the  six  equations  (15)  -(16) 

after  using  (22),  (33)  and  (57)  contain  only  six  unknown  quantities  Fp  Fj,  e^  and  k|  when 

the  externally  applied  forces  and  moments  are  prescribed.  Solution  of  this  set  of  equations 

enables  one  to  use  kj  in  the  curvature-orientation  angle  relations  such  as  (31)  to  determine 

the  latter  i.e.  0i.  Now  lij  (^)  are  known  and  one  proceeds  to  determine  (pi  from  (21)  and 
Xi  from  (9)  using  the  appropnate  boundary  conditions.  The  solution  for  the  first  iteration 

is  complete.  One  enters  (52)  with  ei  and  kj  and  computes  as  and  the  iteration  proceeds. 
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